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We present an analytical solution of the delocalization transition that is induced by an imaginary 
vector potential in a disordered chain [N. Hatano and D. R. Nelson, Phys. Rev. Lett. 77, 570 (1996)]. 
We compute the relation between the real and imaginary parts of the energy in the thermodynamic 
limit, as well as finite-size effects. The results are in good agreement with numerical simulations 
for weak disorder (mean free path large compared to the wavelength). 
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In a recent Letter,EJ Hatano and Nelson have demon- 
strated the existence of a mobility edge in a disor- 
dered ring with an imaginary vector potential. A non- 
Hermitian Hamiltonian containing an imaginary vector 
potential arises from the study of the pinning of voig 
tices by columnar defects in a superconducting cylinder .□ 
Their discovery of a delocalization transition in one 
and twOj^imensional systems has generated considerable 
interest J3ia since all states are localized by disorder in 
one and two dimensions if the vector potential is real. 
Localization in this specific kind of non-Hermitian quan-. 
tum mechanics is referred to as "directed localization" ,0 
because the imaginary vector potential breaks the sym- 
metry between left-moving and right- moving particles, 
without breaking time-reversal symmetry. 

The analytical results of Ref. [l] consist of an expres- 
sion for the mobility edge plus a solution of the one- 
dimensional problem with a single impurity. Here we 
go further, by solving the many-impurity case in one di- 
mension. Most of the technical results which we will need 
were derived previously in connection with the problem 
of localization in the presence of an imaginary scalar po- 
tential. Physically, these two problems are entirely differ- 
ent: an imaginary vector potential singles out a direction 
in space, while an imaginary scalar potential singles out 
a direction in time: A negative imaginary part of the 
scalar potential corresponds to absorption and a positive 
imaginary part to amplification. One might surmise that 
amplification could cause a delocalization transition, but 
in fact all states remain localized in one dijacnsion in the 
presence of an imaginary scalar potential.EnZI 

Following Ref. |l| we consider a disordered chain with 
single-particle Hamiltonian 

(1) 

The operators cj and c^- are creation and annihilation 
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operators, a is the lattice constant, and w the hopping 
parameter. The random potential Vj is chosen indepen- 
dently for each site, from a distribution with zero mean 
and variance . For weak disorder (mean free path much 
larger than the wavelength) , higher moments of the dis- 
tribution of Vj are not relevant. The Hamiltonian is non- 
Hermitian because of the real parameter h, corresponding 
to the imaginary vector potential. The chain of length 
L is closed into a ring, and the problem is to determine 
the eigenvalues e of Ti. If e is an eigenvalue of TL, than 
also e* is one — because H. is real. Real e correspond to 
localizad states, while complex e correspond to extended 
states.EJ 

To solve this problem, we reformulate it in terms of the 
2x2 transfer matrix M^(e) of the chain, which relates 
wave amplitudes at both ends.u The energy e is an eigen- 
value of Ti if and only if Mh(s) has an eigenvalue 1. The 
use of the transfer matrix is advantageous, because the 
effect of the imaginary vector potential is just to multiply 
M with a scalar, 

Mh{e) = e''^Afo(e). (2) 

The energy spectrum is therefore determined by 

det [1 - e^^Moie)] = 0. (3) 

Time-reversal symmetry implies det Afp = ^L. Hence the 
determinantal equation (^ is equivalent tollil 

tr Afo(e) = 2cosh/iL. (4) 

We seek the solution in the limit L ^ oo. 

Since Mq is the transfer matrix in the absence of 
the imaginary vector potential {h = 0), we can use 
the results in the literature on localization in conven- 
tional one-diraensional systems (having an Hermitian 
Hamiltonian) .□ The four matrix elements of A/q are given 
in terms of the reflection amplitudes r, r' and the trans- 
mission amplitude t by 

(A//o)ii = -(lA)det5, (A//o)i2 = r'A, 

(A//o)2i = -r/t, (A/o)22 = l/t, ^''> 

where det S = rr' —t^ is the determinant of the scattering 
matrix. (There is only a single transmission amplitude 



1 



because of time-reversal symmetry, so that transmission 
from left to right is equivalent to transmission from right 
to left.) The transmission probability T — |ip decays 
exponentially in the large-L limit, with decay length ^: 



lim L-^\nT = C 



(6) 



The energy dependence of ^ is known for weak disorder, 
such that |fc|^ ^ 1, where the complex wavenumber k is 
related to e by the dispersion relation 



£ —w cos kc 



(7) 



For real k, the decay length is the localization length ^q, 
given byt2l 



^0 = a{w/uf' sin^(Re ka). 



(8) 



(Since is of the order of the mean free path the con- 
dition of weak disorder requires £ large compared to the 
wavelength.) For complex fc, the decay length is shgrter 
than regardless of the sign of Imfc, according toEl 



2|Imfc| 



(9) 



We use these results to simpHfy Eq. (Q). Upon taking 
the logarithm of both sides of Eq. (^, dividing by L and 
taking the limit L ^ oo, one finds 

\h\~\C^ = Mm L-^\n\l-detSl (10) 

where we have used that L^^lnf L~^ln|/| as i ^ 
oo for any complex function f{L). For complex fc, the 
absolute value of det S is either < 1 (for Im A: > 0) or > 1 
(for Imfc < 0). As a consequence, ln|l — det5| remains 
bounded for L ^ oo, so that the right-hand-side of Eq. 
(I^) vanishes. Substituting Eq. (||), we find that complex 
wavenumbers k satisfy 



\lrak\ = \h\ - l^^K 



(11) 



Together with the expression (||) for the localization 
length this is a relation between the real and imagi- 
nary parts of the wavenumber. Using the dispersion re- 
lation (^, and noticing that the condition |fc|^ 3> 1 for 
weak disorder implies |ImA;| <^ |Refc|, we can transform 
Eq. (|l^) into a relation between the real and imaginary 
parts of the energy, 



llmel 



\h\a^/w^ 



(Re£)2 



2^w2 - (Ree) 



(12) 



The support of the density of states in the complex 
plane consists of the closjsd curve (|l^) plus two line seg- 
ments on the real axis,t3 extending from the band edge 
±w to the mobility edge ±£c- The real eigenvalues are 
identical to the eigenvalues at /i = 0, up to exponentially 
small corrections. The energy Ec is obtained byjputting 
Im£ = in Eq. (p^), or equivalently be equating^ 2^o to 
l/\h\, hence 
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FIG. 1. (a) Data points: eigenvalues of the Hamiltonian 
(0), for parameter values ha — 0.1, u/w = 0.3, and for five 
values of L/a. Dashed curves: analytical large-L limit, given 
by Eq. (^. (Except for the case L — 4000a, spectra are offset 
vertically and only eigenvalues with Ime > are shown.) (b) 
Variance of the imaginary part of the eigenvalues as a function 
of the sample length, for Ree « and for the same parameter 
values as in (a). The data points are the numerical results for 
1000 samples. The solid line is the analytical result (|l5|). 
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(13) 



The delocalization transition at Sc exists for \h\ > he = 

In Fig. the analytical theory is compared with a 
numerical diagonalization of the Hamiltonian (|l|). The 
numerical finite-L results are consistent with the large-L 
limit (dashed curve). To leading order in fluctua- 
tions of Im e around the large-L limit (|l^) are governed 
by fluctuations of the transmission probability T. (Fluc- 
tuations of L~^lnT are of order while the other 
fluctuating contributions to Eq. (^) are p£ order .) 
The variance of In T for large L is knownjij 

2 T 

var InT^ — + 8L|Imfc|e'^«''l^'"'=lEi(-4eo|Imfc|), (14) 
where Ei is the exponential integral. Equating |Im k 
thus 



\h\ + iL-MnT, we find var|ImA:| = i^-^var InT and 



var llmel = 



[l + 27e2''Ei (-27)] [v? - [Reef] , 

(15) 



2L^o 



where 7 = 2|/i|^o — 1- In Fig- |l|b we see that Eq. (|1E 
agrees well with the results of the numerical diagonaliza- 
tion. The fluctuations Alme are correlated over ra, range 
ARee which is large compared to Alme itself, their 
ratio Alme/ARee decreasing oc L^^^^. This explains 
why the complex eigenvalues for a specific sample ap- 
pear to lie on a smooth curve (see Fig. 0a). This curve 
is sample-specific and fluctuates around the large-i limit 

®- 

In conclusion, we have presented an analytical theory 
for the delocalization transition in a single-channel disor- 
dered wire with an imaginary vector potential. We find 
good agreement with numerical diagonalizations, both 
for the relation between the real and imaginary parts of 
the energy in the infinite-length limit and for the finite- 
size effects. In the absence of the imaginary vector po- 
tential, the transfer matrix approach used in this paper 
has been very successful for the study of localization in 
disordered wires with more than one scattering channel. 
We expect that such an extension of the theory is possible 
for non-Hermitian systems as well. 

Discussions with K. B. Efetov and D. R. Nelson moti- 
vated us to do this work. We have also benefitted from 
discussions with T. Sh. Misirpashaev. Support by the 
Dutch Science Foundation NWO/FOM is gratefully ac- 
knowledged. 

Note added: Just before submitting our paper, we 
learned of a different analytical approach to this prob- 
lem by Janik et al.Ej 
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